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Abstract. We calculate weighted mosaic correlations between the recently published Planck 
2015 foreground maps - both anomalous microwave emission (AME) maps, free-free emission, 
synchrotron radiation and thermal dust emission. The weighting coefficients are constructed 
taking account of the signal-to-error ratio given by the data product. Positive correlation is 
found for AME compared with thermal dust emission as well as synchrotron radiation. We 
find AME and free-free emission tending to be anti-correlated, however, when investigating 
different scales, their relationship appears to be more complex. We argue that dust particles 
responsible for AME are pushed out of hot zones in the interstellar medium (ISM). 
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1 Introduction 

Following the release of the PLANCK full mission data, cosmology and astro-particle physics 
are confronted with major challenges concerning the verification of the theories of inflation[l- 
4]. The rapid expansion at the earliest times of the Universe produces gravitational waves. 
This cosmological gravity wave background can be probed by precise measurements of the 
B-mode of polarization of the Cosmic Microwave Background (CMB), as done by [5, 6]. 
Also future effort will be related to this mission by a variety of space- and ground-based, 
experiments (both ongoing and planned) i.e. BICEP3/KECK, SPIDER, EBEX etc. which 
will be able to detect B-mode polarization with unprecedented pixel sensitivity. 

All these experiments have to seriously take into account the recent Planck result [7] 
that there are no zones in the sky, including the BICEP2/KECK region, which are free of 
contamination from galactic origin. In addition, Anomalous Microwave Emission (AME) - a 
component not well understood and only recently added to foreground separation algorithms 
- further hampers one in obtaining a clean CMB map in both, intensity and polarization. 

In order to prepare for coming data, a number of fundamental and practical investiga¬ 
tions have to be performed. The development of new strategies and methods for the removal 
of foregrounds - thermal dust emission at frequencies v > 100 GHz as well as synchrotron 
radiation, free-free emission and AME at v < 100 GHz - is urgently required. 

In this Paper we would like to analyze the morphology of the AME template in relation 
to other foreground components from the Planck 2015 data release [8]. In order to clarify the 
physics of AME, we cross-correlate its template with those of free-free emission, synchrotron 
radiation and thermal dust emission. For this we will use the mosaic correlation method 
(MCM), discussed in [9-11] with some modifications, described below. The specialty of the 
MCM is that it allows one to investigate correlations on different scales. In our analysis we 
restrict ourselves to regions outside the galactic plane in order to place focus on the non-trivial 
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correlations. Furthermore, we weight the data determined by the signal-to-noise ratios of the 
respective maps. It is not surprising that we detect very strong correlations of AME with 
thermal dust and synchrotron emission, widely discussed in the literature [12-14]. However, 
we find significant anti-correlations between AME and free-free emission, especially on small 
scales 1°). This can shed new light on the physics of AME. For larger scales 7°) the 
results are strongly dependent on the weights chosen. 

The structure of the paper is the following: Section 2 will briefly describe the production 
of the foreground templates by the Planck team. In Section 3 we will introduce the coefficient 
of correlation used in this analysis and as well as the weights used to account for errors. In 
section 4 we present the results of our analysis, with and without weighting the data. Lastly, 
section 5 discusses the outcome of our analysis. 

2 The Component Maps 

The Planck 2015 data release provided us with a selection of foreground maps in temperature 
and polarization of which we select five temperature maps for our analysis. In order to provide 
the reader with details about the physical processes behind respective components, we recap 
the essential assumptions for the calculations of the templates by the Planck group in [8]. 1 

Synchrotron radiation is emitted by charged particles of relativistic velocities (mostly 
electrons from cosmic rays) accelerated in the interstellar, galactic, magnetic field. The spec¬ 
trum arising from this radiation is assumed to have unbroken power law form for frequencies 
above 20 GHz. The Planck team used the GALPROP code [16-18] to estimate the spec¬ 
trum caused by cosmic ray electrons. A shift parameter v p acting as a free parameter for the 
subsequent optimization process was introduced to allow for the use of GALPROP at low 
frequencies. While this resulted in the determination of a template function / s (z/), the 408 
MHz Haslam map [19] was used to determine the amplitude A s of the brightness temperature 
for synchrotron radiation: 



( 2 . 1 ) 


where z/, z/q and a are the frequency, the reference frequency for the Haslam 408 MHz map, 
and the shift parameter (a = 0.26) respectively. 

AME, a phenomenon detected by [13-15], currently is assumed to be sufficiently well 
described by spinning dust - rotating dust grains with an electric dipole moment which emit 
radiation at microwave frequencies, as proposed by [20-22]. This model [23-25] was used 
by the Planck team to fit the data. Effectively altering the width of the peak in the AME 
spectrum, the template was chosen to contain two spinning dust components, rather than 
only one, in order to fit the data properly: One had the freedom to change peak frequency 
over the entire sky only constrained by a Gaussian prior centered around 19 GHz with a 
width of 3 GHz, the other kept its peak frequency constant over the sky but could choose it 
freely fitting the data the best which resulted in a peak frequency of 33.35 GHz. The two 
components were evaluated at 22.8 GHz and 41.0 GHz respectively, providing the two AME 
maps (from now on referred to as AME1 and AME2). It was further noted that large sys¬ 
tematic uncertainties correlated with synchrotron and free-free emission are associated with 

^ote, that not all formulas here are presented in their standard form but are rather adopted from [8]. 
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the spinning dust component fit. 


Free-free emission arises from hot electrons being accelerated in intergalactic gas and 
regions of ionized hydrogen. The corresponding antenna temperature Sff (in /jlKrj) for 
free-free emission is given by 

Sff (0 = lo 6 T e (1 - e T// ), Tff = 5.468 • 1(T 2 T~ z ' 2 vf • EM ■ g ff (T e ), (2.2) 


9ff ( T e ) = log 




T 4 = T/IO^K, ug = u/lGHz 


where Tfj is the optical depth for free-free emission, gff is the corresponding Gaunt factor and 
EM is the emissivity measure. The model chosen for the template used by Planck includes two 
free parameters, the emissivity measure EM given by the line-of-sight integration f LOS dl 
and the electron temperature T e . In Planck’s component separation algorithm Commander 
priors were chosen as log(EM) ~ Uniform(—oo, + 00 ) and T e ~ N(7000 K ± 500 K). 


The thermal dust model used is a one component gray body model with a varying 
spectral index: 


sdiy) = A d 


f 1+l3d exp (jug) ~ 1 

\Vd) exp(7i/)-l’ 


h 

kTd 


(2.3) 


where Ud = 545 GHz, (3d is the spectral index, T'd is the dust temperature, Ad is an amplitude 
(including terms of i.e. the black body spectrum as well as the optical depth), h is the Planck 
constant and k is the Boltzman constant. The priors were chosen as (3d ~ iV(1.55 d= 0.1) and 
Td ~ N(23 =b 3K). It is noted that the cosmic infrared background will be absorbed by this 
model as well. 


All discussed maps can be seen in figure 1. Unlike the others, the free-free emission map 
comes in units of pc cm ~ 6 for the emissivity measure than in K for the brightness temperature. 
However, as can be seen from eq. 2.2, since Tff is very small, Sff(v) is approximately 
proportional to tjj, and thereby proportional to EM. Therefore, the morphology of the 
temperature map is hardly any different from that of the EM map which justifies the use of 
the EM map for our correlations. 

All foreground maps come with various extensions. Amongst them we find the map 
obtained by a maximum likelihood method within the Commander algorithm, a mean value 
map and a RMS map including the standard deviations for the mean value map. In our 
analysis, we will use the maximum likelihood map for calculations not including the errors. 
When regarding the error map, however, for obvious reasons, we will use the mean value map. 


3 Mosaic Correlation 

In general, a correlation of two maps in pixel space measures common spatial properties. 
However, these properties might only show up on certain scales. Recently, in a work by 
Hensley and Draine [26] correlations between Planck’s AME and parameters of thermal dust 
emission were analyzed for the masked sky, based on a maximum likelihood approach. The 
conclusion given in this paper, that AME most likely originates from magnetic dust rather 
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Figure 1 . The Planck 2015 foreground templates. From left to right and top to bottom: AME1, 
AME2, synchrotron radiation, thermal dust emission and free-free emission. Here, the maximum 
likelihood maps are shown. 


than spinning dust, cannot be ignored in future work regarding foreground estimation. The 
MCM, as mentioned in sec. 1, provides a means of mapping out regions of correlation on 
different scales. Since we also attempt to weight the data with appropriate coefficients, related 
to the errors given by the data product, we introduce the weighted correlation coefficient. 
Further, we also describe a method of assessing the significance of found correlations suiting 
our purposes. 

3.1 Weighted Mosaic Correlation 

The Mosaic Correlation method computes the correlation between two maps in defined re¬ 
gions, which collectively cover the entire sky without overlapping. The Healpix [27] pixeliza- 
tion offers a convenient way of doing so: pixels at lower resolution than the original resolution 
of the maps mark those regions in which the correlation should be computed. Only one value 
of the correlation coefficient is obtained per "mother pixel". Thereby a full map of correlation 
values is created. Choosing a different resolution of course changes the size of the mother 
pixels and thus the scale at which the correlation maps are investigated. In the following we 
denote a mother pixel by IT Eq. 3.1 defines the weighted and unbiased sample covariance 
between two signals Si and S 2 in IT Wi(p) stands for the weight of signal i in pixel p. 

cov(q,si,s 2 ) = m) ■ v^i (p)«>2(p) (Aw -s^m) (Aw m 

= _ Epgq Vw!(p)w2(p) _ 

(ZpgQ V W 1 (p)w 2 (p)) - EpgQ Wl(p)w 2 (p) 

f(ft) reduces to its well known form of f(fl) = 1/ (Nq — 1) if we set wi{p) = 1, Vp G 12, Vi, 
where Nq is the number of pixels within ft. With Si iW (Q) we denote the weighted average of 
Si(p) in 12, defined as: 


E pe q Wi 


(3.2) 
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Finally, the weighted and unbiased correlation coefficient K in D can be calculated as: 



cov(n, Si, S 2 ) 


(3.3) 


Note, that due to the normalization, the monopole component of the respective maps 
has no impact on K - only the structure in the sky is evaluated. 

3.2 Weights 

Comparing the provided RMS maps with the respective signals reveals that the errors are of 
considerable order. Especially the free-free component’s relative uncertainty is very high at 
high galactic latitudes. We adopt for our analysis a weighting coefficient whose rigor can be 
tuned by choice of the parameter n: 


1 


(3.4) 



c7i(p ) denotes the RMS value for signal Si(p). Very high values of cri(p)/Si(p) will result in a 
weight close to 0 whereas low ratios approach 1. Furthermore, larger values of n correspond 
to weights chosen with greater rigor. A very high n will effectively produce a mask only 
regarding regions of a signal-to-error ratio larger than 1. 

In our following analysis of weighted correlations we will restrict ourselves to n — 1,2 
and as an exemplary choice of a high n, n — 10. Figure 2 shows maps of the resulting weights 
for each of the choices of n for every foreground. We notice that free-free emission has by 
far the lowest signal-to-error ratio (disregarding the galactic plane). The weights of the two 
AME signals show slight domination of the errors at the ecliptic poles, whereas the thermal 
dust and the synchrotron map appear to be the cleanest (Note the ranges of the color scale of 
the respective maps). For another view on the distribution of weights for the different choices 
of n, see figure 3. 

3.3 Assessment of correlations. 

To assess the significance of found correlations we define a statistical ensemble of realiza¬ 
tions, based on general properties of the maps Si and S 2 . We would like to emphasize that 
both these maps represent statistically inhomogeneous and non-Gaussian signals. Hereto, all 
foreground’s PDFs can be seen in the left panel of figure 4. 2 

In this paper, we focus on the properties of the AME maps in comparison to the syn¬ 
chrotron, free-free and thermal dust emission maps. We therefore create artificial signals 
based on the respective AME maps, as described below, which we then correlate with the 
remaining foreground maps for comparison with the real correlation. Firstly, we decompose 
the AME signal into spherical harmonics: 


max 


S ame {6,<)>) = Y, E ISimle^YUe^) 


(3.5) 


1=0 m=—l 


2 In order to make the PDFs of all foregrounds well visible in one plot, we actually show the PDFs of 
the foreground’s temperature relative to its median value. Additionally, we exclude the 10% pixels with the 
highest values as well as the 1% pixels with the lowest values. 
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Figure 2. Maps of the weights calculated for the respective foregrounds. From top to bottom: AME, 
AME2, free-free emission, synchrotron emission and thermal dust emission. From left to right: n = 1, 
n = 2 and n = 10 



Figure 3. PDFs of the values of the weights for the respective foregrounds. From left to right: n = 1, 
n = 2 and n = 10. The legend for all plots can be found in the middle panel. 
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where \Si m \ are the amplitudes and <F/ m are the phases of each (Z,m) mode, respectively. The 
power spectrum of S AME -map is then given by C AME — ^ |S) 7m | 2 , independent 

of the phases <F/ m . For a statistically inhomogeneous and non-Gaussian signal the statistical 
distribution of phases is far from a uniform distribution, with different correlations between 
(Z,m) and (Z + /c,m + n) modes for various k,n = 0,1, 2... . At the same time, detecting a 
common morphology in two maps, say AME and another foreground map, means that the 
phases 4>/ m and the phases for the foreground signal, T/ m , are strongly correlated, ~ 

The last statement gives us an idea of how to generate a statistical ensemble of realizations of 
artificial AME maps, which could only have chance correlations with synchrotron, free-free or 
thermal dust maps. For that, we extract the amplitudes of the AME maps \Si m \ from Eq(3.5) 
and select the phases randomly from a uniform distribution within [0, 2 tt] . Combining 
IS/mI and 1> lm to SZ nd = |5, m | exp(i^ m ) for each realization. The right panel in figure 4 
shows the PDFs of all 1000 random AME1 maps, generated as just described, compared to 
the actual signal. 3 The same process has been applied to the AME2 map (not shown here). 




Figure 4. Left: Probability density functions for used foreground templates. Right: Probability 
density function of AME1 compared with that of a randomly generated AME1 sky as described in 
the text. 


4 K-maps for AME and synchrotron, free-free and dust emission maps 

In this section we summarize the results of mosaic cross-correlations for the two AME skies 
with synchrotron, free-free and thermal dust maps. N S id e = 256 for all input maps which de¬ 
fines the area of each pixel to q p = Att/N Pi where N p is the number of pixels for N^e • The sizes 
of the mother pixels, D, are chosen according to the lower resolutions N' side = 64, 32, 16, 8 
which correspond to the sizes Qq — N$iq p , with Nq — 16, 64, 256, 1024. Nq is the number 
of pixels within one mother pixel for a given N' side . The angular extent of each pixel in the 
original resolution is 9 p ~ 14 arcmin and for the mother pixels 0' p ~ 55 arcmin, 1.8°, 3.7° and 
7.4°. 

Before showing the results of the weighted correlations, we choose to neglect the errors 
of both maps to obtain a regular, unweighted correlation. For this, we use the maximum 
likelihood map rather than the mean value map, as already mentioned in sec. 2. In the 

3 Again for better visibility, we show PDFs of the temperature T relative to the RMS of T. Since the random 
realizations have both, negative and positive values, for this plot, we consider the absolute value of the maps. 
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subsequent section we then include the errors and weight the mean value data with the three 
different coefficients introduced in section 3.2, corresponding to different levels of rigor towards 
the data. 

4.1 Unweighted correlations 

We start our analysis with unweighted correlations, corresponding to the assumption that 
(Ji — 0 and thus — 1. Figure 5 shows the correlation maps for the combinations of AME1 
and AME2 with free-free emission, synchrotron radiation and thermal dust emission for all 
discussed pixelizations. Strong correlations are found for all combinations along the galactic 
plane. When showing the PDFs of each of the maps in figure 6, we therefore plot both - a 
PDF of all pixels of the map as a dotted line and another as a solid line after excluding the 
galactic plane, simply taken to be a uniform band of 30° opening angle centered around the 
galactic plane. We notice that, as expected, high correlation values are removed. In addition, 
the average PDF over 1000 simulations, computed as described in section 3.3, are shown as 
a black line with those of all single simulations in gray indicating the spread. 



Figure 5. From top to bottom: Unweighted correlation maps between AME1 and free-free emission, 
AME2 and free-free emission, AME1 and synchrotron radiation, AME2 and synchrotron radiation, 
AME1 and thermal dust emission, AME2 and thermal dust emission. From left to right: Q contains 
1024, 256, 64 and 16 pixel. Here, the maximum likelihood maps are used. The cqs were all assumed 
to be zero. 







Figure 6. PDFs of the unweighted correlation coefficients K for different combinations of foregrounds 
and each of the Q sizes. From top to bottom: Correlation between AME1 and free-free emission, 
AME2 and free-free emission, AME1 and synchrotron radiation, AME2 and synchrotron radiation, 
AME1 and thermal dust emission, AME2 and thermal dust emission. Colored lines give distributions 
of K’s from the actual foregrounds, solid for masked sky, dotted for unmasked sky. The gray lines 
show all distributions of K’s resulting from random simulations of AME1 or AME2 correlated with 
the respective foregrounds. Black lines are averages over all 1000 random realizations. Simulations 
are shown for the masked sky. From left to right: Q contains 1024, 256, 64 and 16 pixels. 
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4.2 Weighted correlations 

Now, we include the errors <Ji(p) and therefore use the mean value maps as Si(p) instead of 
those of maximum likelihood. We saw in figure 3.4 that the weights corresponding to the 
free-free map will have by far the largest impact on the correlations, compared to the other 
foregrounds’ weights. Indeed, we see that after taking the weights under consideration only 
the correlations of AME with free-free emission deviate significantly from the unweighted 
ones. Therefore, we here only show the results of those correlations, while we move the 
results regarding synchrotron radiation and thermal dust emission to the appendix (figures 
14 - 19). The correlation maps for n = 1,2 and 10 can be seen in figures 7, 8 and 9, while 
the corresponding PDFs are shown in figures 10, 11 and 12, respectively. 



Figure 7. From top to bottom: Weighted-Correlation maps between AME1 and free-free emission, 
AME2 and free-free emission. From left to right: Q contains 1024, 256, 64 and 16 pixel. Here, we 
weighted with Wi(n = 1). 



Figure 8. Same as figure 7, but here, weighted with weighting coefficients Wi(n = 2). 






Figure 9. Same as figure 7, but here, weighted with weighting coefficients Wi(n = 10). 
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Figure 10. PDFs of weighted correlation coefficients K (with Wi(n = 1)) for different combinations 
of foregrounds and each of the D sizes. From top to bottom: Correlation between AME1 and Free- 
free emission, AME2 and Free-free emission. Colored lines give distributions of K’s from the actual 
foregrounds, solid for masked sky, dotted for unmasked sky. The gray lines show all distributions 
of Ks resulting from random simulations of AME1 or AME2 correlated with the free-free emission 
signal, and their averages (black line). Simulations are shown for the masked sky. From left to right: 
D contains 1024, 256, 64 and 16 pixels. 



Figure 11. Same as figure 10, but here, weighted with weighting coefficients Wi(n = 2) 






Figure 12. Same as figure 10, but here, weighted with weighting coefficients Wi(n = 10) 
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4.3 Comparison of weighted and unweighted results 

It is apparent that the simulations as well as the actual cross-correlations change their dis¬ 
tribution function according to the size of the mother pixels D. In order to portray this fact 
and also to compare all correlations, in figure 13 we simply show the mean value of the co¬ 
efficients E[K ) = K • P(iF), where P(K) is the PDF for a given resolution. E(K) acts 
as the average similarity of the AME components in respect to the other foregrounds from 
the smallest angular scale of about 0 P ~ 0.9° up to 9 P « 7.4°. The corresponding range of 
multipoles for those scales can be estimated as l • 6 P < 1, which gives l < 126 for Nq = 16 
and l < 16 for Nq = 1024. 

As one could already see in the PDFs (figures 6,10-12,17-17) and now more clearly from 
figure 13, the correlations of foreground components can be separated into correlations, more 
pronounced at small (as AME and thermal dust) and at large angular scales (as AME and 
free-free). The most distinct positive correlations ( E(K ) ^ 10 — 50%) are detected for AME1 
and AME2 with thermal dust at all scales, increasing for larger scales. Hardly influenced by 
scale, the correlations between AME and synchrotron emission he steady around 10 —20 %. As 
already mentioned, for both, AME with thermal dust and AME with synchrotron, weighting 
does not change their behavior noticeably. 

The unweighted correlation between AME and free-free emission shows an opposite trend 
to that of the thermal dust correlation. Here E(K ) ranges between -20 and -40% for small 
scales and goes to zero for larger ones. When applying weights to the data, the strong negative 
correlation is pushed to higher values. For large scales the values even grow to be positive up 
to +20%. 

5 Conclusion 

We have computed weighted mosaic correlations between recently released foreground tem¬ 
plates of the Planck Collaboration. This method allows for the determination of correlations 
within defined regions, filling out the entire sky. Since the method only takes place in pixel 
space, it can be applied to inhomogeneous, anisotropic signals on the sky, as is the case here. 
When altering the size of the pixels, one obtains correlations corresponding to different scales. 
We then simply chose the expectation value of the distribution of correlation coefficients for 
each probed scale as a rough estimator of the manifestation of correlation and thereby ob¬ 
tained a scale dependency of the correlation. Furthermore, we included the errors on the 
signals, provided by the Planck data product, by introducing appropriate weights. These 
could be tuned to vary the emphasis put on data points with signal-to-error ratios below or 
above 1. 

A strong positive correlation was detected for AME and thermal dust, for both, the 
unweighted and weighted case. Also the correlation of both AME components with syn¬ 
chrotron radiation give a pronounced positive correlation. As mentioned before, the positive 
correlation with thermal dust is in accordance with the earliest studies of AME. Due to the 
common relation of both - the dust particles causing AME and electrons causing synchrotron 
radiation - to the galactic magnetic fields, the correlation with synchrotron emission was also 
to be expected. 

Both AME components show a strong anti-correlation with free-free emission at small 
angular scales for the unweighted case, especially at the smallest angular scales. We noticed 
an even stronger anti-correlation for AME2 with free-free emission than for AME1. After 
including the weights in the calculation, the anti-correlation decreases noticeably on small 
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Figure 13. The mean value of the cross-correlation coefficients E(K) for all the combinations of 
foregrounds. Top left: unweighted correlations, top right: weighted correlations with Wi(n = 1), 
bottom left: weighted correlations with Wi(n = 2) and bottom right: weighted correlations with 
Wi(n = 10). The legend for all plots can be found in the upper right panel. Dashed lines correspond 
to values from the entire sky, whereas solid lines correspond to the masked sky. 


scales, while we even see positive values on larger scales. Negative correlation is restored again 
on all scales for the choice of the most rigorous weights at n = 10. The weighted correlations 
differ substantially from the unweighted ones when we regard the free-free template. Due 
to its very small signal-to-error ratios on medium to high galactic latitudes, the weights 
Wi(n = 1,2) seem to induce positive correlations among AME and free-free emission, formerly 
anti-correlated. Values of those pixels with a low signal-to-error ratio are pushed down towards 
zero by the weights. We therefore obtain positive correlations with the AME skies because 
they as well have low values in exactly those regions. The choice of Wi(n — 10) basically 
forces all pixels with signal-to-error ratios less than 1 completely to zero. It appears that 
those pixels causing the positive correlation for the choice of n = 1,2 are now contributing 
to the negative correlation. Our conclusion must be that the positive correlations at high 
galactic latitudes should not be taken seriously. Moreover, we should acknowledge that the 
negative correlations closer to the galactic plane are robust as they survive even the most 
aggressive weights. 

We interpret these results by claiming that the nano-particles causing AME are pushed 
out of the hot zones within our galaxy. The electrons, responsible for free-free emission, 
remain in these zones, perhaps being held back by abundant ions. 
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Understanding the composition and the dynamics of the ISM is essential in separating 
galactic foregrounds and in the long run is vital for a correct insight in the early Universe, 
including the detecting of a cosmological gravitational wave background. Assuming that 
the AME skies provided by the Planck Collaboration describe the true anomalous emission 
sufficiently well, we believe that our result should be taken into account in future work on 
models of AME and the ISM. 

Acknowledgment s 

We are very appreciative of helpful comments by K. Gorski and A. Frejsel and thank espe¬ 
cially P. Naselsky for his support and useful discussions in this work. 

This work is based on observations obtained with Planck 4 , an ESA science mission with in¬ 
struments and contributions directly funded by ESA Member States, NASA, and Canada. 
The development of Planck has been supported by: ESA; CNES and CNRS / INSU-IN2P3- 
INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); 
CSIC, MICINN and JA (Spain); Tekes, AoF and CSC (Finland); DLR and MPG (Germany); 
CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ire¬ 
land); FCT/MCTES (Portugal); and PRACE (EU). A description of the Planck Collaboration 
and a list of its members, including the technical or scientific activities in which they have 
been involved, can be found at the Planck web page. 5 

We acknowledge use of the HEALPix 6 package [27]. This work is supported by Danmarks 
Grundforskningsfond which allowed the establishment of the Danish Discovery Center, the 
Villum Foundation in support of the Deep Space project, the National Natural Science Foun¬ 
dation of China (Grant No. 11033003), the National Natural Science Foundation for Young 
Scientists of China (Grant No. 11203024) and the Youth Innovation Promotion Association. 


4 http://www.esa.int/Planck 

5 http://www. cosmos, esa.int/web/planck/planck-collaboration 

6 http://healpix.sourceforge.net 


- 14 - 




A Additional Figures 



Figure 14. From top to bottom: Weighted-Correlation maps between AME1 and synchrotron radi¬ 
ation, AME2 and synchrotron radiation, AME1 and thermal dust emission, AME2 and thermal dust 
emission. From left to right: Q contains 1024, 256, 64 and 16 pixel. Here, we weighted with Wi{n = 1). 



Figure 15. Same as figure 14, but here, weighted with weighting coefficients Wi(n = 2). 
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Figure 17. PDFs of weighted correlation coefficients K (with Wi{n = 1)) for different combinations of 
foregrounds and each of the Q sizes. From top to bottom: AME1 and Synchrotron radiation, AME2 
and Synchrotron radiation, AME1 and Thermal Dust emission, AME2 and Thermal Dust emission. 
Colored lines give distributions of K’s from the actual foregrounds, solid for masked sky, dotted for 
unmasked sky. The gray lines show all distributions of K’s resulting from random simulations of 
AME1 or AME2 correlated with the respective foregrounds. Black lines are averages over all 1000 
random realizations. Simulations are shown for the masked sky. From left to right: Q contains 1024, 
256, 64 and 16 pixels. 
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Figure 18. Same as figure 17, but here, weighted with weighting coefficients Wi{n = 2) 
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Figure 19. Same as figure 17, but here, weighted with weighting coefficients Wi(n = 10) 

References 

[1] Guth, A. H. The Inflationary Universe: A Possible Solution to the Horizon and Flatness 
Problems , Phys.Rev. D , 23 (1981) 347-356 

[2] Mukhanov, V. F. and Chibisov, G. V., Quantum Fluctuation and Nonsingular Universe (In 
Russian), JETP Lett 33 (1981) 532-535 

[3] Linde, A. D., A new Inflationary Universe Scenario: A possible Solution of the Horizon, 
Flatness, Homogeneity, Isotropy and Primordial Monopole Problems , Phys.Lett. B , 108B 
(1982) 389 

[4] Albrecht, A. and Steinhardt, P. Cosmology for Grand Unified Theories with Radiatively 
Induced Symmetry Breaking , PRL , 48 (1982) 1220 

[5] Bicep2 Keck Array Collaboration, Ade, P. et ah, Detection of B-Mode Polarization at Degree 
Angular Scales by BICEP2 , Phys. Rev. Lett. 112 (2014) 241101, arxiv: 1403.3985 

[6] BICEP2/Keck Array and Planck Collaborations, A Joint Analysis of BICEP2/Keck Array and 
Planck Data , arxiv: 1502.00612 

[7] Planck Collaboration Int. XXX, Planck intermediate results. XXX. The angular power spectrum 
of polarized dust emission at intermediate and high Galactic latitudes , A&A (2014) in press 


- 19 - 






































































[8] Planck Collaboration X, Planck 2015 results. X. Diffuse component separation: Foreground 
maps , arxiv: 1502.01588 

[9] Hansen, M. et al., Faraday Rotation as a diagnostic of Galactic foreground contamination of 
CMB maps , MNRAS 426 (1) (2012) 57-69, arxiv:1202.1711 

[10] Verkhodanov, O.V. and Khabibullina, M.L., Dominant Multipoles in WMAP5 Mosaic Data 
Correlation Maps , Bull. Spec. Astrophys. Obs ., 65 (2010) 390 

[11] Verkhodanov, O. V., Khabibullina, M. L., Majorova, E. K. Tessellated mapping of cosmic 
background radiation correlations and source distributions , Astrop. Bull ., 64, 3 (2009) 263-269 

[12] Kogut, A., Banday, A. J., Bennett, C. J., Gorski, K. M., Hinshaw, G. and Reach, W. T. 

High-Latitude Galactic Emission in the COBE Differential Microwave Radiometer 2 year sky 
maps , ApJ, 460, (1996) 1-9 

[13] Leitch, E. M., Readhead, A. C. S., Pearson, T. J. and Myers, S. T., An anomalous component 
of galactic emission , ApJ 486 (1997) L23-L26 

[14] de-Oliveira-Costa, A., Kogut, A., Devlin, M. J., Netterfield, C. B., Page, L. A. and Wollack, E. 
J. Galactic Microwave Emission at Degree Angular Scales , ApJL , 482, (1997) L17-L20 

[15] Kogut, A. et ah, A mechanism of non-thermal radio-noise origin , ApJL 464 (1997) L5 

[16] Moskalenko, I. V. and Strong, A. W., Production and Propagation of Cosmic-Ray Positrons 
and Electrons , ApJ , 493, (1998) 694 

[17] Strong, A. W., Moskalenko, I. V. and Ptuskin, V. S., Cosmic-Ray Propagation and Interactions 
in the Galaxy , Annual Review of Nuclear and Particle Science , 57, (2007) 285 

[18] Orlando, E. and Strong, A. W., Galactic synchrotron emission with cosmic ray propagation 
models , MNRAS , 436 (2013) 2127 

[19] Haslam, C. G. T. et ah, A 408 MHz all-sky continuum survey. II - The atlas of contour maps , 
Astron. Astrophys. Suppl. Ser. 47 (1982) 1-142 

[20] Erickson, W. C. A mechanism of non-thermal radio-noise origin , ApJ 126 (1957) 480 

[21] Draine, B. T. and Lazarian A., Diffuse Galactic Emission from Spinning Dust Grains , ApJL 
494 (1998) L19 

[22] Draine, B. T. and Lazarian A., Electric Dipole Radiation from Spinning Dust Grains , ApJ 508 
(1998) 157 

[23] Ali-Haimoud, Y., Hirata, C. M. and Dickinson, C.A refined model for spinning dust radiation , 
MNRAS , 395, 2 (2009) 1055-1078 

[24] Silsbee, K., Ali-Haimoud, Y. and Hirata, C. M .Spinning dust emission: the effect of rotation 
around a non-principal axis , MNRAS , 411, 4 (2011) 2750-2769 

[25] Ali-Haimoud, Y. 2010, SpDust/SpDust.2: Code to Calculate Spinning Dust Spectra, 
Astrophysics Source Code Library 

[26] Hensley, B. S. and Draine, B. T, A case against spinning PAHs as the source of the Anomalous 
Microwave Emission , arxiv: 1505.02157 

[27] Gorski, K., Hivon, E., Banday et ah, HEALPix: A Framework for High-Resolution 
Discretization and Fast Analysis of Data Distributed on the Sphere , ApJ 622 (2005) 759 


- 20 - 



